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PREFACE 


The  work  reported  herein  was  conducted  by  the  Arnold  Engineering 
Development  Center  (AEDC),  Air  Force  Systems  Command  (AFSC). 

The  work  was  done  by  ARO,  Inc.  {a  subsidiary  of  Sverdrup  &  Parcel 
and  Associates,  Inc. ),  contract  operator  of  AEDC,  AFSC,  Arnold  Air 
Force  Station,  Tennessee.  The  work  was  done  under  ARO  Project 
No.  VD227,  and  the  manuscript  (ARO  Control  No.  ARO-VKF-TR-73 - 
135)  was  submitted  for  publication  on  October  1,  1973. 
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1.0  INTRODUCTION 


1.1  ,  PROBLEM  DEFINITION 

Brayton  cycle  refrigeration  systems  are  frequently  used  in  plant 
installations  because  of  component. simplicity  and  acceptable  efficiency. 
Such  refrigeration  systems  normally  consist  of  either  a  centrifugal  or  a 
piston  compressor,  heat  exchangers,  and  an  expander.  The  most 
familiar  application  of  the  Brayton  cycle  is  the  gas  turbine.  The  ideal 
processes  of  this  cycle  are  isentropic  compression,  isobaric  heat  trans¬ 
fer,  and  isentropic  expansion  followed  by  isobaric  heat  transfer.  The 
ideal  processes  of  the  Brayton  cycle  refrigerator  are  usually  isothermal 
compression,  isobaric  heat  transfer,  and  isentropic  expansion.  Despite 
the  difference  in  processes,  most  literature  refers  to  the  refrigeration 
cycle  described  either  as  a  modified  Brayton  cycle  or  simply  as  a 
Brayton  cycle. 

Refrigerators  installed  in  plants  are  normally  designed  for  efficient 
steady-state  operation.  Brayton  cycle  refrigerators  are  also  frequently 
designed  for  flight  use,  but  they  have  some  problems  not  dealt  with  in 
this  study.  Refrigerators  designed  for  efficient,  steady-state  operation 
and  a  minimum  first  cost  of  installation  may  not  possess  optimum  cool¬ 
down  characteristics.  When  a  refrigerator  is  used  to  provide  the  heat 
sink  for  a  large  aerospace  testing  facility,  this  additional  time  may  be 
very  expensive.  As  an  example,  a  typical  test  run  of  12  hrs  may  require 
7  hrs  of  refrigerator  operation  to  cool  the  chamber  to  the  required  tem¬ 
perature. 

An  analytical  and  experimental  study  was  implemented  with  the  goal 
of  minimizing  the  time  required  to  cool  a  large  panel  in  a  vacuum  cham¬ 
ber  to  20°K. 

1.2  SOLUTION  APPROACH 

The  approach  taken  in  minimizing  cooldown  time  was  both  analytical 
and  experimental.  A  preliminary  analytical  approach  was  established 
from  first  principles.  As  the  cooldown  of  an  existing  system  was  moni¬ 
tored,  modifications  and  additions  to  the  analytical  model  were  made. 

The  study  consisted  of  (1)  developing  an  optimization  function,  related 
to  cooldown  time,  which  could  be  used  to  compare  performance, 

(2)  writing  a  computer  program  to  calculate  the  optimization  function, 

(3)  utilizing  the  program  to  perform  a  sensitivity  analysis  of  opera¬ 
tional  and  hardware  variables,  (4)  running  comparative  tests  with  an 
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operational  refrigerator  to  test  operational  variables  against  analytical 
predictions,  and  (5)  evaluating  the  effects  of  possible  cycle  variations. 

2.0  ANALYTICAL  STUDY 
2.1  THE  OPTIMIZATION  FUNCTION 

A  single  term  was  desired  to  compare  the  results  of  a  wide  variety 
of  modifications  to  the  basic  cycle  variables.  The  time  required  to 
cool  the  load  from  ambient  temperature  to  some  arbitrary  value  in  the 
test  region  was  the  obvious  choice.  However,  an  exact  calculation  of 
time  for  a  particular  load  can  become  more  involved  than  is  justified 
for  this  analysis.  Instead,  an  optimization,  inversely  proportional  to 
time,  was  developed.  This  function  was  used  to  compare  the  relative 
effects  of  operational  and  design  parameters  for  a  fixed  load  mass  and 
material. 

Figure  1  is  a  schematic  of  the  refrigeration  system,  and  Fig.  2 
shows  a  temperature -entropy  diagram  of  the  system  at  a  particular 
operating  point.  The  numbers  listed  on  both  figures  will  be  used  as 
subscripts  in  the  equations  to  be  derived.  These  subscripts  will  denote 
particular  locations  in  the  thermodynamic  cycle  (hereafter  referred  to 
simply  as  cycle). 


Figure  1.  Refrigeration  system  schematic. 
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Figure  2.  Temperature-entropy  diagram. 


To  derive  the  optimization  function,  a  control  volume  of  gas  was 
selected  inside  the  load  (Fig,  3).  The  thermal  energy  added  to  the  gas 
as  it  passes  through  the  control  volume  over  an  increment  of  time  is 
expressed  as 

t 

Q  =  /m  Cp(T  12  -  Tn)dt 
o 

where 

Q  =  thermal  energy 
t  -  time 

m  =  mass  flow  rate 

Cp  =  specific  heat  at  constant  pressure 
T  -  temperature 
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This  same  amount  of  thermal  energy  is  taken  from  the  mass  of  the 
load  (the  metal  tube  through  which  the  gas  passes),  lowering  its  temper¬ 
ature  as  follows: 

,Tf 

Q  =  /  WCdT 

T. 

1 


Both  the  mass  flow  rate  and  the  gas  temperatures  Tn  and  T12  are 
unknown  functions  of  time.  They  are  known  functions  of  temperatures 
and  are  set  by  cycle  performance,  as  will  be  seen  later.  The  specific 
heat  of  helium,  the  flow  medium  considered  in  this  study,  does  not  vary 
significantly  between  300  and  20°K  at  the  pressures  considered  (less 
than  150  psi).  For  a  parti cular  load  mass,  W,  the  thermal  energy,  Q, 
to  be  removed  can  be  calculated  over  any  desired  temperature  incre¬ 
ment,  since  the  metal  specific  heat  variation  with  temperature  is  known. 


If  an  average  mass  flow,  m,  and  gas  temperature  difference 
(T12  -  T11)  is  used  over  a  short  increment  of  time,  the  two  equations 
shown  previously  can  be  equated  and  solved  as  follows: 

t  Tf 

/ m  Cp(Ti2  -  Tu)dt  =  /  WCdT 
o  T. 


m(Ti2  -  Tn)Cp  /  dt 


W  f  CdT 

T. 

1 


At  =  W  /  CdT/rh(Ti2  -  Tn)Cp 

T. 

1 
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For  a  given  system,  the  integral  f  CdT  is  a  fixed  number.  Thus, 

T 

1 

to  minimize  the  time  required  to  cool  the  system  over  the  increment 
where  in  and  (T^2  ”  Tn)  are  appropriately  constant,  it  is  only  neces¬ 
sary  to  maximize  the  product  of  (m)  x  (T^2  ~  Tn).  This  product  must 
be  taken  over  many  small  temperature  increments  during  the  cooldown 
process  to  be  valid,  since  both  m  and  T12  •  Tn  vary  considerably.  It 
is  necessary  to  form  a  numerical  integral  of  all  these  products  through 
the  cooldown  process  to  determine  a  measure  of  the  heat  removed  by 
the  gas.  Cycle  variations  that  maximize  this  integral  also  cause  the 
cooldown  time  to  be  minimized.  The  integral 


T12f 

/  m(T12-Tn)dT 

T12i 


is  designated  as  the  optimization  function,  OPF.  To  assure  accuracy, 
T12  was  stepped  down  in  increments  of  only  2°K.  The  corresponding 
value  of  Tn  was  determined  from  cycle  calculations  and  is  described 
in  Section  2.  2  of  the  present  report. 

It  should  be  noted  that  the  derivation  just  described  does  not  con¬ 
sider  any  external  heat  leak  into  the  load.  It  was  assumed  that  all  the 
refrigeration  went  to  cooling  the  load.  Typical  heat  leaks  are  calcu¬ 
lated  for  an  existing  system  to  justify  this  assumption.  The  load  con¬ 
sisted  of  approximately  1200  ft^  of  aluminum  panel  weighing  approxi¬ 
mately  4500  pounds.  The  panel  was  located  inside  a  chamber  with 
pressures  less  than  10”^  torr.  Other  panels  cooled  to  80°K  surrounded 
the  load.  All  panels  were  painted  black  (emissivities  assumed  =  1.0). 
These  surrounding  panels  are  normally  at  operating  temperature  when 
the  refrigerator  is  started.  Thus,  when  a  cooldown  is  started,  thermal 
energy  is  lost  by  radiation.  As  the  temperature  of  the  load  is  cooled 
below  80°K,  the  load  begins  to  receive  thermal  energy  by  radiation 
from  the  surrounding  panels.  Figure  4  is  a  plot  of  the  radiation  rate 
of  the  load  as  a  function  of  panel  temperature.  For  comparison,  a 
typical  refrigeration  curve  as  a  function  of  the  temperature  of  the 
panel  was  plotted  in  Fig.  5.  Radiation  is  seen  to  be  mostly  insignifi¬ 
cant  to  the  refrigeration  available. 
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Figure  4.  Load  radiation. 


Figure  5.  Refrigerator  capacity. 


Conduction  heat  transfer  into  the  load  can  easily  be  made  insignifi¬ 
cant  by  proper  design.  With  a  tensile  stress  of  10,  000  psi,  a  stainless 
steel  rod  with  a  cross  section  of  0.45  in.2  can  support  the  4500-lb 
panel.  Assuming  the  rod  is  only  1  ft  long  and  connects  solidly  to  a 
300°K  wall,  use  of  the  proper  thermal  conductivity  integral  shows  the 
heat  leak  to  be  only  2.  89  w. 

2.2  CYCLE  ANALYSIS 


The  value  of  ~  is  set  by  the  performance  capability  of  the 
cycle  components.  In  the  refrigeration  system  studied  in  this  analysis 
the  compression  process  was  carried  out  with  a  three -stage  piston 
compressor  with  water  intercooling.  Each  stage  was  equipped  with 
water-jacketed  cylinder  walls  and  was  represented  in  the  analysis  as  a 
polytropic  process.  Heat  exchange  was  essentially  an  isobaric  process, 
and  the  expander  was  corrected  from  the  ideal  isentropic  process  by  an 
expander  efficiency. 
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Equations  for  the  various  parameters  were  derived,  using  T12  as 
an  input  variable  to  ultimately  calculate  Tji  and  thus  T^g  ■  Tn.  A 
value  for  m  was  input  as  a  function  of 

With  the  expander  bypassed,  only  the  liquid  nitrogen/gaseous 
helium  (LN2/GHe)  heat  exchanger  provides  cooling.  A  parameter  de¬ 
scribing  the  performance  of  such  a  heat  exchanger  is  the  heat  exchanger 
effectiveness.  This  parameter  is  defined  as  the  actual  heat  transferred 
divided  by  the  maximum  possible  heat  transfer.  For  the  L^/GHe  heat 
exchanger  the  effectiveness  is 


13  14 

LN2  "  T13  -  TLN2 

Provisions  were  made  to  input  heat  leaks  as  well  as  pressure  drops 
into  the  cycle.  Heat  leaks  were  entered  as  temperature  changes  for  con¬ 
venient  insertion  into  the  cycle  calculations.  The  modified  effectiveness 
equation  for  the  L^/GHe  heat  exchanger,  considering  heat  leak,  is 

eLN2  =  <t13  -  T14  +  AT13-14>/(T13  ”  TLN2) 

Similarly,  the  effectiveness  of  the  GHe/GHe  heat  exchanger,  with  heat 
leak,  is 


eGHe  =  (T16  '  T15  +  ATi5-16)/(T7  -  T15) 

Using  the  above  equations,  the  sequence  of  calculations  for  T^j 
with  the  expander  bypassed  is 

t13  =  t12  +  AT12-13 

t14  =  T13{1  -  eLN2)  +  tLN2  eLN2  +  AT13-14 
t15  =  T14  +  t14-15 

T8  =  €GHeT15  +  T7(x  -  eGHe>  +  AT7-8 
T7  =  300°K  by  water  cooling 

t16  *  U  -  eGHe^T15  +  T7  eGHe  “  AT15-16 
Tg  =  T8  +  AT8_g 

T10  =  T9  +  at9-10 
T11  =  T10  +  ATio-11 
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When  the  expander  is  operating  and  the  LN2/GHe  heat  exchanger  is 
being  utilized,  all  temperatures  except  T^q  are  calculated  just  as  before. 
The  temperature  leaving  the  expander,  Tjq,  can  be  calculated  from  the 
pressure  ratio  across  the  expander  and  the  expander  efficiency,  as 
follows  (note  that  k  =  Specific  heat  ratio): 

Tio  -  T9{l  -  €exp  [l  -  <P9/Pi0)^]}  +  AT9_10 

where 

p9  =  p7  “  (AP7-8  +  AP8-£)) 

P7  =  Pi(P2/Pl)3 

p10  =  P1  +  ap16-1  +  aP13-1G  +  aP14-15  +  AP13-14 
+  AP12-13  +  AP11-12  +  AP10-11 
When  the  LN2/GHe  heat  exchanger  was  not  being  utilized. 


t14  =  t13  +  at13-14 

A  computer  program  was  written  for  performing  the  calculations 
listed  above  on  a  digital  computer  using  a  FORTRAN  IV  language.  In¬ 
put  data  were  the  system  parameters  such  as  stage  pressure  ratio, 
P2P1,  heat  leaks,  AT,  pressure  drops,  AP,  and  various  efficiencies. 
Also,  codes  were  established  for  utilizing  the  proper  equations  de¬ 
pending  on  the  particular  configuration.  This  was  necessary  to  account 
for  operational  variables  such  as  TeXp,  the  temperature  of  the  gas  com¬ 
ing  from  the  load  (T^)  at  which  the  expander  is  added  to  the  cycle  (°K), 
and  Tstop,  the  temperature  of  the  gas  coming  from  the  load  (T^  at 
which  the  LJN^/GHe  heat  exchanger  is  removed  from  the  cycle  (°K).  It 
was  assumed  in  this  study  that  the  L^/GHe  heat  exchanger  would 
always  be  in  operation  at  the  start  of  a  cooldown.  The  variable  Tst0p 
was  used  so  that  the  computer  would  use  the  proper  equations  after  the 
LN2/GHe  exchanger  was  removed  from  the  cycle.  It  is  always  re¬ 
moved  from  the  cycle  at  some  point  because  it  begins  to  add  thermal 
energy  to  the  refrigerant  rather  than  removing  it.  The  expander 
could  be  brought  on  at  any  temperature,  TeXp,  including  ambient,  but 
it  was  assumed  that  it  would  remain  operational  for  the  rest  of  the  cool¬ 
down. 

During  the  cooldown  process  T12  was  reduced  2°K  at  a  time.  The 
mass  flow  was  a  program  input  as  a  function  of  T^.  For  each  value  of 
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T^'2  ;a' value  of  Tn  was  calculated  utilizing  the  appropriate  cycle  equa¬ 
tions.  The  product  m(T^2  “  Tn)  was  then  formed  and  stored  as  a  func¬ 
tion  of  Ti2-  When  T12  was  reduced  from  300  to  20^,  these  values  of 
m(Ti2  -  T11)  were  taken  from  storage,  and  the  optimization  function, 

T12f 

OPF  =  f  m(T  12  -  Tn)dT 
T12± 

was  calculated  numerically,  using  the  trapezoidal  rule. 

2.3  RESULTS  OF  THE  COMPUTER  ANALYSIS 

To  obtain  some  basic  knowledge  about  the  system,  it  was  desirable 
first  to  determine  the  parameters  that  had  the  most  significant  effects 
on  cooldown  time.  The  more  important  parameters  could  then  be 
studied  in  more  detail  with  the  goal  of  improving  performance.. 

Examination  of  the  cycle  equations  revealed  nine  parameters  per 
study,  as  follows:  (1)  compressor  stage  pressure  ratio;  (2)  L^/GHe 
heat  exchanger  effectiveness;  (3)  GHe/GHe  heat  exchanger  effective¬ 
ness;  (4)  expander  efficiency;  (5)  pressure  drop  upstream  of  the  ex¬ 
pander;  (6)  pressure  drop  downstream  of  the  expander;  (7)  the  temper¬ 
ature  at  which  the  expander  is  turned  on;  (8)  the  temperature  at  which 
the  LN2/GHe  heat  exchanger  is  turned  off;  and  (9)  the  temperature  of 
the  LN2.  Since  an  operational  refrigerator  was  available  for  compari¬ 
son  with  the  analytical  model,  reasonable  starting  values  for  each  param¬ 
eter  were  available.  These  basic  values  were  as  follows: 

1.  P2P1  =  2.  04 

2.  eLN2  =  0*  98 

3-  eGHe=0*988 

4.  eexp  =  0.  73 

5.  APn_12  =  2.0  psi 

6.  AP-7_g  =2.0  psi 

7.  Texp  =  100"K 

8.  Tgtop  “  100°K 

9.  TLN2  =  80°K 
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A  systematic  series  of  computer  calculations  was  performed  in 
which  each  parameter  was  varied  above  and  below  the  basic  value  while 
all  other  parameters  were  held  constant  at  the  basic  value.  During  this 
first  set  of  runs  it  was  assumed  that  the  piping  in  the  system  would 
cause  sufficient  pressure  losses  to  cause  a  variation  in  mass  flow  dur¬ 
ing  cooldown.  It  was  known  from  previous  observations  of  the  operating 
system  that  full  mass  flow  was  possible  through  the  load  only  at  tem¬ 
peratures  below  approximately  20°K.  An  equation  was  derived  for  de¬ 
termining  the  mass  flow  at  higher  temperatures  utilizing  the  Panning 
Equation.  This  did  not  give  a  wide  variation  with  mass  flow,  and  the 
resulting  computer  studies  indicated  that  the  expander  should  be  turned 
on  at  300°K  to  minimize  cooldown  time.  Early  experimental  runs 
showed  the  derived  mass  flow  equation  to  be  incorrect,  and  the  cool¬ 
down  was  actually  slowed  by  bringing  the  expander  on  at  300°K.  Study 
of  mass  flow  measurements  to  the  chamber  were  not  meaningful  from 
a  quantitative  standpoint  because  of  unresolved  instrument  error.  How¬ 
ever,  it  did  become  apparent  that  the  flow  was  choked  somewhere  in 
the  expander.  An  expression  was  derived  which  could  be  used  to  calcu¬ 
late  the  mass  flow  through  the  expander  as  a  function  of  cycle  conditions. 
The  equation  for  mass  flow  into  the  expander  is 


rfi  =  pAV 


where  p  is  the  gas  density,  A  is  the  flow  area,  and  V  is  the  flow  ve¬ 
locity.  Utilizing  the  equation  of  state  and  the  equation  for  the  velocity 
of  sound,  the  equation  becomes 

m  =  ^  x  A  x  yj kgRT 

where  R  is  the  gas  constant,  k  is  the  specific  heat  ratio,  and  g  is  the 
gravitational  constant.  Pressures  and  areas  did  not  change  when  the 
expander  was  brought  on  line.  The  only  variable  in  the  above  equation 
is  T  or  Tg,  the  temperature  of  the  gas  entering  the  expander.  Thus, 
the  form  of  the  mass  flow  equation  is 

K 

m  =  /= 

VTg 

with  the  constant  K  evaluated  to  be  1.  565  for  m  =  0.  35  lbm/sec  and 
Tg  =  20°K. 

Utilizing  this  mass  flow  variation,  each  parameter  was  again  sys¬ 
tematically  varied  above  and  below  the  nominal  value.  Results  of  these 
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calculations  are  shown  in  Figs.  6  through  13.  Five  of  the  parameters, 
expander  efficiency,  LN2/GHe  heat  exchanger  effectiveness,  GHe/GHe 
heat  exchanger  effectiveness,  pressure  drop  after  the  expander,  and 
pressure  drop  before  the  expander,  are  seen  to  be  quite  linear.  Com¬ 
pressor  stage  pressure  ratio  is  a  quadratic  function  but  only  varies 
0.04  percent  over  the  reasonably  assumed  range  of  operation.  A 
straight-line  curve  fit  of  OPF  versus  P2P1  was  made  using  least 
squares  and  the  computer  output.  As  might  be  expected  from  closer 
examination  of  the  cycle  equations,  the  effects  on  the  optimization  func¬ 
tion  of  LN2  temperature  variation  and  the  temperature  at  which  the 
L^/GHe  heat  exchanger  is  turned  off  are  coupled.  The  variation  in 
optimization  function  with  varying  expander-on  temperatures  is  signifi¬ 
cantly  nonlinear  but  is  conveniently  represented  by  a  quadratic  function. 


Figure  6.  Effect  of  varying 

expander  efficiency. 


LNji'GHe  Heat  Exchanger  Effectiveness 

Figure  7.  Effect  of  varying  LN2/GHe  heat 
exchanger  effectiveness. 
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Figure  8.  Effect  of  varying  GHe/GHe  heat 
exchanger  effectiveness. 
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Figure  9.  Effect  of  varying  pressure 
drop  after  expander. 
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Figure  10.  Effect  of  varying  pressure 
drop  before  expander. 


Figure  11.  Effect  of  varying  compressor 
stage  pressure  ratio. 
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Tstop  Expander-On  Temperature,  °K 

Figure  12.  Effect  of  varying  Tstop  and  Figure  13.  Effect  of  varying  expander-on 
LIM2  temperature.  temperature. 

2.4  THE  OBJECTIVE  FUNCTION 

To  study  the  effects  each  variable  has  on  cooldown  time,  it  is  de¬ 
sirable  to  have  a  common  basis  of  comparison.  To  utilize  the  method¬ 
ology  developed  in  the  field  of  operations  research  for  optimization  of 
the  system's  performance,  a  single  equation  describing  the  effects  of 
each  variable  is  required.  This  single  equation  that  includes  all  the 
variables  is  called  an  objective  function.  In  this  study  the  objective 
function  becomes  an  equation  that  is  used  to  calculate  the  optimization 
function  as  a  function  of  the  hardware  and  operational  variables.  The 
same  symbol,  OPF,  was  used  for  the  objective  function.  In  the  follow¬ 
ing  explanation,  xj  for  i  =  1  to  n  is  used  as  a  substitute  for  the  variables 

eLN2’  etc*  ’  ^or  brevity- 

For  any  selected  set  of  hardware  and  operational  variables  the  ob¬ 
jective  function  is 

OPF  =  OPF(Basic)  +  d(OPF) 
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where  d(OPF)  can  be  written 


d(OPF)  = 


9(OPF) 

9x^ 


dx^ 


9(OPF) 

3x2 


dx2  + 


9(OPF) 
9x 


9(OPF) 

dxi 


is  recognized  as  the  slope  of  the  plots  of  OPF  versus  vari¬ 


able  i  found  in  Figs.  6  through  13.  The  amount  a  particular  variable  is 
changed  from  Basic  is  dx^.  The  amount  the  objective  function  is  changed 
from  Basic  because  of  all  variable  changes  is  d(OPF). 


With  the  variables  set  at  their  basic  values  (Section  2.  2), 
OPF(Basic)  is  8292.3001. 


As  an  example,  two  variables  will  be  used  to  demonstrate  how  the 
objective  function  was  derived. 


OPF 

d(OPF) 


=  OPF(Basic)  +  d(OPF) 
90PF  J  90PF 

'^exp  eXp  ^LN2 


deLN2 


d(OPF)  =  649.  2(eexp-0.  73)  +  8101(eLN2  -  0.  98) 

OPF  =  8292.  3001  +  649.  2  eexp  +  8101  eL]vj2  -  8412 
OPF  =  649.  2  eexp  +  8101  €ln2  “  120.  6 


which  agrees  within  0.  2  percent  of  the  values  calculated  by  the  computer 
program  for  eexp  =  0.  80  and  ^LN2  =  0.  99. 

The  technique  described  above  was  used  with  linear  variations  of 
OPF  versus  x^.  For  nonlinear  variations  such  as  OPF  versus  Tg^p,  a 
different  curve-fitting  technique  was  used.  Figure  13  was  represented 
by  the  following  equation,  holding  all  other  variables  constant: 

OPF  =  8269  -  0.  036  T2 

exp 


As  an  example,  this  equation  will  be  combined  with  the  equation 
previously  derived  for  varying  eexp  and  £LN2’ 

OPF  =  649.  2  epxn  +  8101  eL N_  -  0.  036  T2  +  k 
exp  -^^2  exp 
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where  k  is  a  constant  evaluated  by  substituting  basic  values  of  OPF, 
eexp»  eLN2J  Texp.  By  using  these  techniques,  the  objective  equa¬ 
tion  was  finally  determined  to  be  as  follows: 

OPF  =  649.  2  eexp  +  8101  eLN2  +  17,  700  e  GHe 

-9.02  APn_12-1.28  AP7_8 

+  231.  8  P2P1  -  0.  069  T2,  +  11.02  T 

stop  stop 

-  68  Tt  -  0.  036  T2  -  12,  585 
J-,1N2  exp 


2.5  OPTIMIZATION 


With  an  objective  function,  the  second  step  in  utilizing  operations 
research  methodology  to  obtain  an  optimum  solution  is  to  select  the 
appropriate  constraints.  In  this  study  the  maximum  stage  pressure 
(P2P1)  is  2.  04  and  is  set  by  hardware  limitations.  Values  for  heat  ex¬ 
changer  effectiveness  (both  eG j.je  and  eLN2^  anc*  the  expander  efficiency 

(eeXp)  must  be  less  than  1.  0  for  thermodynamic  reasons.  All  pressure 
drops  obviously  must  be  greater  than  zero.  With  the  refrigeration 
process  starting  at  ambient  conditions,  the  expander  must  be  turned  on 
and  the  LN2/GHe  heat  exchanger  turned  off  at  some  temperature  less 
than  300°K.  The  liquid  nitrogen  temperature  under  normal  conditions 
will  be  greater  than  77°K  since  that  is  its  boiling  temperature  at  ambi¬ 
ent  pressure.  However,  the  system  used  in  the  experimental  study 
was  pressurized,  and  80°K  was  used  as  the  basic  temperature.  With 
the  reasons  given  above,  the  constraints  are  listed  below: 


P2P1  <  2.04 
eGHe  -  1-° 
eLN2  -  1*° 


eexp  — 
AP7-8  - 


Texp 

^stop 

TLN2 


< 

< 

< 


1.0 

0 

300°K 

300°K 

80°K 
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Some  of  the  above  conditions  could  have  been  predicted  without 
such  optimization  procedures.  However,  intuition  does  not  tell  the 
amount  of  influence  each  variable  has  on  the  objective  function  when 
conditions  vary  from  optimum.  Determining  these  quantities  is  re¬ 
ferred  to  as  a  sensitivity  analysis.  The  variation  in  the  objective 
function  with  a  change  in  variable  i  is  3(OPF)/9xj.  For  each  variable 
considered,  the  sensitivity  value  is  listed  below. 


3(opf) 

a(P2Pi) 

=  231.8 

3(0  PF) 
9<eGHe> 

=  17, 700 

3  (OPF) 
9<eLN2) 

=  8101 

3  (OPF) 

9  (eexp) 

=  649.2 

3  (OPF) 
3(AP7_8) 

=  -1.  28 

3  (OPF) 

=  -9.02 

9(AP11-12) 

3(OPF) 

9(Texp) 

=  -0,072  TeXp 

3(OPF) 

9(Tstop) 

=  -0.  069  Tstop+  11.02 

3(OPF) 

9<tLN2> 

=  -68 

The  relative  effects  of  these  variables  can  be  expressed  in  more 
meaningful  terms  by  recalling  that  OPF  was  derived  to  be  inversely 
proportional  to  time.  For  the  system  that  was  tested,  a  typical  cool¬ 
down  took  approximately  7  hrs  with  OPF  =  8292  with  eQHe  =  0.  988. 
Using  the  sensitivity  coefficient,  the  predicted  change  in  OPF  with  a 
drop  of  ^QHe  *°  9.  987  would  be 

AOPF  =  17,  700(0.  988  -  0.  987)  =  17.  7 
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The  predicted  cooldown  time  would  be 


8292-17. 7 
8292  X 

A  corresponding  change  of  0.01 


hrs  =  6.  98  hrs 
in  £eXp*  in  time,  is 


AOPF  =  649.  2(0.  98  -  0.  97)  =  6.  49 
8292  -  6.49 


time  = 


8292 


x  7  =  6.  99  hrs. 


3.0  TEST  RESULTS 

During  the  conduction  of  a  series  of  tests  in  the  Mark  I  Aerospace 
Chamber  at  AEDC,  cooldown  data  were  taken  for  comparison  with  ana¬ 
lytical  studies.  The  Mark  I  Chamber  is  an  84-ft-high  by  42-ft-diam 
vacuum  vessel  entirely  lined  with  panels  which  are  cooled  to  77°K. 
Pressures  are  approximately  10”°  torr  during  testing.  A  100-ft-long 
by  6-ft-wide  aluminum  panel  is  suspended  inside  the  liner  and  cooled 
to  less  than  20°K  with  gaseous  helium.  Refrigeration  is  supplied  with 
a  Brayton  cycle  refrigerator,  and  helium  from  the  refrigerator  flows 
through  passages  in  the  aluminum  panel.  Cycle  variables  are  approxi¬ 
mately  those  listed  previously  as  Basic. 

Platinum  temperature  sensors  were  attached  to  the  panel  on  the 
inlet  and  exhaust  manifolds.  During  the  cooldown  different  procedures 
were  used  to  check  assumptions  made  during  analysis.  A  typical  plot 
of  cooldown  data  is  shown  in  Fig.  14.  It  should  be  noted  that  it  is  ex¬ 
tremely  difficult  to  obtain  what  might  be  called  laboratory-type  data 
on  such  a  large  system  while  it  is  undergoing  tests  for  some  other  pur¬ 
pose.  There  was  a  separate  circuit  which  was  also  cooled,  but  it  took 
only  a  small  amount  of  refrigerant.  It  was  turned  on  at  various  un¬ 
known  times  during  the  cooldown  process.  Heat  leaks  were  found  dur¬ 
ing  some  cooldowns  and  were  eliminated.  One  heat  leak  was  actually 
located  by  referring  to  computer  printouts  of  predicted  temperatures 
and  checking  against  some  installed  temperature  indicators.  No  indi¬ 
cators  were  installed  at  the  heat  leak  point.  These  problems  did  not 
prevent  the  cooldown  from  being  completed  as  required,  but  it  did 
make  comparisons  with  analytical  work  less  accurate. 

Cooldown  began  with  the  2500-lb  aluminum  GHe/GHe  heat  ex¬ 
changers  at  various  temperatures.  These  heat  exchangers  are  very 
well  insulated,  and  if  only  a  few  days  pass  between  cooldowns,  they 
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stay  cold.  The  usual  operational  procedure  is  to  run  the  refrigerator 
in  a  bypass  mode  around  the  chamber  to  cool  the  entire  refrigerator, 
including  the  heat  exchangers,  prior  to  starting  cooldown  on  the  load. 
However,  this  is  not  always  possible.  In  the  computer  program  which 
calculates  the  optimization  function,  the  warm  heat  exchanger  is  simu¬ 
lated  as  a  heat  leak. 


Figure  14.  Typical  load  cooldown. 

Table  I  is  a  tabulation  of  the  operating  conditions  and  cooldown 
times  observed  during  tests.  Only  data  for  cooldowns  which  ran 
smoothly  are  shown.  Measured  temperatures  of  73  and  78°K  for  Ts^0p 
are  probably  from  5  to  7°K  lower  than  the  actual  gas  temperature  but 
were  used  as  the  point  to  take  the  L^/GHe  heat  exchanger  off  the  line. 

Table  1 .  Comparison  Between  Analysis  and  Test 
Conditions  Test  Analysis 


^exp» 

°K 

Tstop* 

°K 

Initial  Heat 
Exchanger 
Temp,  °K 

Percent 
Change  in 

Time  for 

Basic  =  6.  9  Hr 

OPF 

Percent 
Change  in 
OPF  from 
Basic 

220 

90 

169 

11.9 

7018 

18 

122 

107 

183 

2.  98 

8181 

1.35 

130 

73 

179 

1.49 

8156 

1.67 

125 

78 

178 

2.8 

8205 

1.  06  . 
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Values  of  OPF  were  calculated  for  the  particular  operating  condi¬ 
tions  using  Eq.  (1),  from  Section  2.  4. 

The  results  presented  confirm  the  predicted  effects  of  changes  in 
the  operational  variables  TeXp  and  Tst0p.  Starting  the  expander  at  a 
high  temperature  slows  the  cooldown  time  because  of  insufficient  mass 
flow.  Varying  Tst0p  over  the  range  from  107  to  73°K  only  varies  the 
overall  cooldown  time  a  small  amount. 

4.0  CYCLE  VARIATIONS 

4.1  DESCRIPTION  OF  CONSIDERED  VARIATIONS 

There  are  many  possible  cycle  variations  that  might  be  considered 
if  this  were  a  study  of  a  proposed  system  rather  than  an  existing  one. 

As  it  is,  only  changes  that  either  overcame  serious  existing  problems 
or  could  be  added  with  a  minimum  of  hardware  modification  are  con¬ 
sidered. 

Both  experimental  test  and  analysis  confirmed  the  limiting  compon¬ 
ent  on  the  test  refrigerator  to  be  the  expander  and  its  inability  to  pass 
full  compressor  mass  flow  at  higher  temperatures.  The  results  of  a 
study  of  the  system  assuming  full  compressor  mass  flow  throughout  the 
cooldown  process  are  given  in  Section  4.  2.  This  could  be  accomplished 
by  adding  an  expander  in  parallel  with  the  existing  one.  After  the  system 
was  cooled  to  such  a  temperature  that  one  expander  could  handle  full 
mass  flow,  one  expander  could  be  shut  down. 

Some  literature  on  cryogenic  refrigerators  assumes  an  additional 
L^/GHe  heat  exchanger  located  in  between  two  GHe/GHe  heat  exchangers. 
Such  an  addition  could  be  provided  for  the  existing  system  and  is  con¬ 
sidered  in  Section  4.3. 

Reheat  cycles  are  frequently  used  to  improve  efficiency.  This 
modification  is  checked  for  transient  performance  in  Section  4.4. 

4.2  ADDING  AN  EXPANDER 

Since  adding  an  expander  would  overcome  serious  difficulties  in  the 
existing  system,  a  complete  analysis  was  run  in  that  configuration.  An 
objective  function  was  derived  using  methods  described  previously, 
assuming  full  compressor  mass  flow  during  the  entire  cooldown.  The 
computer  output  is  plotted  in  Figs.  15  through  22.  The  resulting  equa¬ 
tion  for  OPF  is  as  follows: 
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OPF  =  1488P2P1  -  8.  141  AP?_8  -  58.  55  APX1.12 
+  13,  179  eGHe  +  4949  eLN2  +  4157  eexp 
+  11-  82  Texp  -  0.  1116  T|top  +  17.  85  Tstop 
-  17,245 


11. 100 


11,000 


10,900 


1G.SG0 


10,700 


10,600 


re  15.  Effect  of  varying  pressure  drop 
before  expander  (m  =  0.35). 


Figure  16.  Effect  of  varying  expander 
efficiency  (m  =  0.35). 


LN2''GHe  Heat  Exchanger  Etlectiveness 

Figure  17.  Effect  of  varying  LN2/GHe 
exchanger  effectiveness 
(m  =  0.35). 


Pressure  Drop  Atler  Expander,  psi 

Figure  18.  Effect  of  varying  pressure  drop 
after  expander  {m  =  0.35). 
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Figure  19.  Effect  of  varying  GHe/GHe  Figure  20.  Effect  of  varying  expander-on 
heat  exchanger  effectiveness  temperature  (m  =  0.35). 

(m  =  0.35) 


Figure  21.  Effect  of  varying  compressor  stage  Figure  22.  Effect  of  varying  TItop 
pressure  ratio  (m  =  0.35).  ( rin  =  0.35). 
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By  comparing  the  sensitivity  terms  (see  Table  2)  for  the  two  differ¬ 
ent  mass  flow  distributions,  one  can  see  that  the  stage  pressure  ratio, 
is  6.4  times  as  sensitive  with  rh  =  0.35  as  it  is  with  m  =  1.  565AJ"T9. 
However,  changes  in  ^LNg  are  approximately  1.3  times  less  sensitive. 

With  this  modification,  cooldown  time  could  be  reduced  24  percent. 


Table  2.  Mass  Flow  Distributions 


Sensitivity  Terms 

m  =  1.  565/  \/Tq 

m  =  0.  35 

aoPF/aP2Pi 

232.4 

1488 

aOPF/aAPbefore  expander 

-1.  279 

-8. 141 

aOPF/3APafter  expander 

-9.  017 

-58. 55 

30PF  /  9eQHe 

17,  770 

13,  180 

aoPF/aeLN2  ' 

6,492 

4,  949 

90PF/  exp 

649 

4,  157 

aOFF/3Texp 

-0.072  Texp 

11.  82 

aoPF/aTstop 

-0.069  Tstop+  11.02 

-0.  223  Tgj-0p  +  !! 

aOPF/3TLN2 

-68 

(Not  considered) 

4.3  ADDING  AN  LN2/GHe  HEAT 

EXCHANGER 

. 

Adding  a  second  LN2/GHe  heat  exchanger  in  the  middle  of  the  exist¬ 
ing  GHe/GHe  heat  exchanger  would  be  a  relatively  simple  task.  In  the 
analytical  model  previously  described,  the  GHe/GHe  heat  exchanger  is 
a  single  unit.  However,  on  the  refrigerator  that  was  tested  there  are 
two  GHe/GHe  heat  exchangers  in  series.  Transient  performance  calcu¬ 
lations  indicated  that  such  a  modification  would  reduce  cooldown  times 
by  only  3.  3  percent. 

4.4  THE  REHEAT  CYCLE 

A  reheat  cycle  involves  using  two  expanders  and  dividing  the  load 
into  two  sections.  The  cold  gas  from  the  first  expander  is  routed  through 
the  first  section  of  the  load  and  back  to  the  second  expander.  Cold  gas 
from  the  second  expander  is  routed  through  the  remaining  part  of  the 
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load.  In  the  transient  analysis  it  was  assumed  that  full  compressor 
mass  flow  would  pass  through  both  expanders  and  that  the  temperature 
of  the  gas  entering  both  expanders  was  the  same. 

This  modification  would  cut  the  cooldown  time  in  half. 

5.0  CONCLUSIONS 

The  conclusions  of  this  study  are  as  follows: 

1.  The  inability  of  the  expander  in  the  GHe  refrigerator  of 
the  Mark  I  Chamber  to  pass  full  compressor  mass  flow 
is  the  system's  greatest  deterrent  to  rapid  cooldown. 

2.  An  additional  expander  will  allow  full  compressor  mass 
flow  to  be  utilized  throughout  the  cooldown  process  and 
should  reduce  cooldown  times  by  24  percent. 

3.  Precooling  the  system's  GHe/GHe  heat  exchanger  prior 
to  test  can  greatly  reduce  cooldown  time. 

4.  Operator  attention  should  be  given  to  maintaining  full 
mass  flow  through  the  load,  particularly  when  the 
expander  is  being  used.  A  more  reliable  gauge  should 
be  installed  to  allow  this  to  be  accomplished. 

5.  With  the  existing  system  the  expander  should  be  utilized 
when  the  gas  temperature  returning  from  the  load 
reaches  115°K. 

6.  With  the  existing  system  the  LN^/GHe  heat  exchanger 
should  be  utilized  until  the  gas  temperature  returning 
from  the  load  reaches  80°K. 

7.  If  future  installations  require  larger  load  masses,  cool¬ 
down  times  may  become  excessive.  A  reheat  cycle 
using  two  expanders  and  existing  equipment  should  re¬ 
duce  cooldown  times  to  half  their  present  value. 
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NOMENCLATURE 


C 


Cp 

g 

k 

m 


OPF 

P 

P2P1 


Q 


R 

T 

T 

exp 


T 


stop 


t 

W 

AP7-8 

AP11-12 


Specific  heat 

Specific  heat  at  constant  pressure 
Gravitational  constant 
Specific  heat  ratio 
Mass  flow  rate 

Optimization  function,  proportional  to  cooldown  time 
Pressure,  psi 

Compressor  stage  pressure  ratio 
Thermal  energy,  or  heat  transferred 
Gas  constant 
Temperature,  °K 

Temperature  of  the  gas  coming  from  the  load  (T^)  at 
which  the  expander  is  added  to  the  cycle, 

Temperature  of  the  gas  coming  from  the  load  (T^)  at 
which  the  LN2/GHe  heat  exchanger  is  removed  from 
the  cycle,  °K 

Time 

Mass 

A  representative  pressure  drop  preceding  the  expander 
A  representative  pressure  drop  following  the  expander 
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eexp 

eGHe 

eLN2 


Expander  effectiveness 

Heat  exchanger  effectiveness,  GHe/GHe  heat  exchanger 
Heat  exchanger  effectiveness,  LN2/GHe  heat  exchanger 
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